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ABSTRACT 


The  integral  equation,  first  obtained  by  H,  Maruo, 
which  determines  the  pressure  distribution  generatina  flow 
past  a  slender  ship  of  vanishing  draft,  is  studied  further. 
New  results  obtained  include  predictions  of  singular  center- 
plane  effects  of  gravity  for  pointed  bodies,  a  similarity 
solution  for  ships  with  cusped  parabolic  waterplanes,  and 
some  preliminary  numerical  solutions  of  the  integral  equation 
in  the  general  case. 


1 .  Introduction 

This  paper  is  an  extension  of  a  study  by  Maruo  (1967)  on  "ships 
of  small  draft",  "flat  ships",  or  "planing  surfaces",  all  of  which 
are  equivalent  descriptions.  The  small-draft  assumption  allows 
linearization  of  the  free-surface  boundary  condition,  as  in  the 
comparable  case  of  thin  ships,  or  ships  of  small  beam  (Michell,  (1898)). 
However  the  present  linearized  problem  is  much  harder  to  solve,  since 
the  generating  singularity  distribution  (effectively  a  distribution  of 
pressure  points  on  the  limiting  waterplane)  is  not  explicitly  given  in 
terms  of  hull  shape,  but  requires  solution  of  an  integral  equation. 

This  problem  is  analogous  to  the  lifting-surface  problem  of  aerodynamics, 
whereas  the  thin-ship  problem  corresponds  to  the  simpler  thickness 
problem  of  aerodynamics. 

Although  we  also  give  here  a  brief  re-consideration  of  the  general 
flat-ship  problem,  to  emphasize  some  aspects  not  discussed  by  Maruo  (1967)  , 
our  attention  in  the  present  paper  is  mainly  devoted  to  the  low-aspect- 
ratio  limit.  Thus  the  wetted  length  of  the  ship  is  supposed  much  greater 
than  its  beam,  the  latter  having  already  been  assumed  much  greater  than 
the  draft  by  the  flat-ship  requirement.  The  ship  is  therefore  not  only 
flat ,  hut  also  slender. 

We  give  an  alternative  derivation  of  an  integral  equation  equivalent 
to  one  obtained  by  Maruo  (1967) ,  which  has  as  its  unknown  function  a 
pressure  distribution  representing  the  ship.  This  integral  equation  is 
also  obtainable  from  the  high-Froude-number  s lender-body  theory  of 
Ogilvie  (1967),  by  assuming  that  the  ship  is  not  only  slender t  but  also 
flat. 

Maruo' s  low-aspect-ratio  flat-ship  integral  equation  is  formally 
valid  only  at  moderately-high  Froude  numbers,  specifically  such  that 
U2B/gL2  is  a  quantity  of  order  unity,  where  U  is  ship  speed,  B  its 
beam  and  L  its  length,  and  g  is  gravity.  The  equation  reduces  to  that 
of  low-aspect-ratio  wing  theory  in  aerodynamics  as  g  +  0  .  One  approach 
to  practical  solution  of  any  planing  problem,  whether  or  not  the  aspect 
ratio  is  low,  is  to  expand  in  an  asymptotic  series  for  Very  high  Froude 
number,  commencing  with  the  aerodynamic  g  =  0  limit  as  the  leading 

( e . sy .  Wang  s  Rispin,  (1971).  Maruo  (1967)  obtains  the  first  two 
terms  in  this  series  for  the  lift  on  a  flat  delta  wing,  and  we  give  here 
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an  alternative  treatment  of  this  class  of  expansion,  for  general  hull 
shapes.  In  particular,  we  demonstrate  very  strong  effects  of  gravity  near 
the  center  plane  of  pointed  bodies. 

We  also  observe  that  at  all  Froude  numbers,  the  low-aspect-ratio 
flat-ship  integral  equation  possesses  a  "similarity"  solution,  such 
that  the  pressure  distribution  has  the  same  shape  at  all  stations.  This 
linearized  but  gravity-dependent  result  should  not  be  confused  with  the 
well-known  conical  similarity  solution  for  non-linear  planing  or  water 
entry  in  the  absence  of  gravity  (Gilbarg,  (I960),  p.  360).  In  fact  the 
present  geometrical  requirement  is  for  a  cusped  parabolic  waterplane  shape 
but  an  arbitrary  section  shape,  whereas  the  non-linear  zero-gravity  solution 
requires  a  triangular  plan  form  and  section  shape. 

The  low-aspect-ratio  flat-ship  integral  equation  is  amenable  to 
direct  computation,  and  we  present  here  some  preliminary  examples  of  its 
numerical  solution.  Much  more  work  needs  to  be  done  to  derive  efficient 
procedures,  and  the  present  computer  program  can  only  l  considered  as  a 
crude  first  attempt.  However,  the  results  are  of  considerable  interest, 
indicating  rather  dramatic  gravity  effects  especially  near  the  center 
plane,  as  predicted  analytically,  and  confirming  Maruo's  (1967)  estimate 
of  the  lift  coefficient  of  a  delta  wing  at  sufficiently  high  Froude  number. 


2.  The  General  Flat  Ship  Problem 

We  use  a  rather  unconventional  co-ordinate  system  (x,y,s),  as  in 
Tuck  and  von  Kerczek  (1968)  and  as  sketched  in  Figure  1.  The  ship  is 
supposed  fixed  with  its  bow  at  s=0  and  stern  at  s=L  in  a  stream  U.  Thus 
the  total  flow  field  velocity  is 

q  =  V(Us  -Hj>)  ,  (2.1) 

where  $  is  the  perturbation  velocity  potential. 

The  body  equation  is 

y  =  Tl(x,s)  ,  (2.2) 

where  H  is  generally  expected  to  be  negative,  "-fl"  being  the  depth 
of  the  buttock  line  x=constant  at  station  s.  Equation  (2.2)  is  supposed 
to  hold  for 

| x |  <  b(s)  ,  (2.3) 

where  b(s)  is  the  half-waterplane  width  at  station  s.  For  |xj  >  b 
we  may  suppose  that  (2.2)  defines  the  water  surface  elevation.  The  hull 
boundary  condition  is 

<J>y  *  (u  +  <}>s)ns  +  <J>xnx  ,  (2.4) 

to  be  applied  on  the  exact  hull  surface  y=n 

We  first  make  the  small-draft  approximation,  introducing  a  small 
parameter  a  measuring  the  draft/length  ratio.  Keeping  only  leading 
order  terms  with  respect  to  a  ,  the  boundary  condition  (2.4)  reduces  to 

<j>y  =  uns  ,  on  y  =  0.  (2.5) 

It  is  important  to  note  that  the  small-  a  approximation  is  a  regular  one, 
as  distinct  from  the  (potentially)  singular  perturbation  represented  by 
the  small-e  slenderness  approximation  to  be  applied  next,  where  e 
measures  the  beam/length  ratio.  We  shall  assume  that  ct<<£  so  that 
(2.5)  may  be  taken  to  hold  quite  accurately  when  we  come  to  make  the  small 
e  approximation . * 


*It  is  of  interest  to  note  that,  according  to  Acosta  and  DeLong  (1971) ,  the 
infinite-Froude-number  slender-planing-surface  analysis  of  Tulin  (1956)  is 
valid  in  the  opposite  limit  £<<a  . 


The  boundary  condition  (2.4)  also  correctly  gives  the  exact  kinematic 
condition  on  the  unknown  free  surface  y=n  for  |x|  >  b  .  This  has  to  be 
supplemented  by  the  dynamic  condition 

l  +  u4>s  +  j  |V4>I2  +  gn  -  o  ,  (2.6) 

if  the  excess  of  pressure  over  atmospheric  at  the  free  surface  is  P;  usually 
P=0.  Again,  the  small-draft  approximation  enables  linearization  not  only 
of  (2.4)  but  also  of  (2.6)  to  give 

2.  +  u<|>s  +  gn  -  0  ,  on  y  -  0,  (2.7) 

which  combines  with  (2.5)  to  give  the  linearized  free-surface  condition 

g<t> y  +  U2<t>ss  =-  -  Ps  .  (2.8) 

p 

If  P=0,  this  reduces  to  the  usual  equation 

g<fy  +  U24>ss  =  0.  (2.9) 

However,  we  shall  generate  solutions  by  means  of  pressure  distributions  P, 
the  velocity  potentials  then  satisfying  (2.8)  whenever  P  ?  0.  Note  that 
(2.9)  results  from  the  small-a  approximation,  and  that  alone;  when  we 
subsequently  take  e  as  small,  (2.9)  may  be  considered  as  exact. 

The  general  flat-ship  problem,  with  E  not  necessarily  small,  is  that 
of  solving  the  full  Laplace  equation 

$XX  +  ^yy  +  ^SS  =  0  (2'10) 

in  the  space  y<0,  subject  to  the  hull  condition  (2.5)  on  the  portion 
|x|<b(s)  of  the  plane  y«0  occupied  by  the  projection  of  the  hull,  and  the 
linearized  free-surface  condition  (2.9)  on  the  portion  |x|>b(s).  In 
addition  we  expect  to  require  some  kind  of  radiation  condition  at  infinity, 
and  a  Kutta-type  condition  that  the  pressure  reduces  to  atmospheric 
pressure  at  any  sharp  trailing  edge  in  order  that  the  free  surface  leave 
such  an  edge  smoothly. 

This  problem  can  be  converted  into  an  integral  equation,  which  is 
the  finite-Froude-number  analogue  of  the  lifting-surface  integral  equation 
of  aerodynamics.  Maruo  (1967)  gives  one  method  for  accomplishing  this; 
perhaps  more  directly  we  may  set  ourselves  the  task  of  finding  an  unknown 
surface  pressure  distribution  P(x,s)  which  generates  the  free-surface 
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displacement  n(x,s).  The  corresponding  integral  connection*  between  P 
and  n  may  be  obtained  from  well-known  formulae,  e.g.  Wehausen  and  Laitone 
1960,  p.  598.  For  example 


(x,s)  = 


-Im  II  d£dap(S,0) 
¥ 


TT/2 

J  d0secO 


(2.11) 


f 


2  -ik(s-O)  cos0 


dk 


k  e 


co s  ( k  ( x-Q  sin0) 


9  9 

k-0-j  secz0 


where  the  path  of  k- integration  goes  above  the  pole. 

We  shall  not  attempt  to  solve  this  integral  equation  here,  since  our 
concern  is  with  the  low-aspect-ratio  case.  However,  several  questions  are 
worth  noting.  Maruo  (1967)  states  that  "the  kernel  of  the  integral  equation 
is  complicated  enough  to  frustrate  any  attempt  at  solving  it."  This  view  is 
perhaps  a  little  too  pessimistic.  The  kernel  is  simply  the  complete  solution 
for  a  travelling  three-dimensional  pressure  point,  and  a  number  of  similar 
computations  have  been  carried  out  on  an  ad  hoc  basis  recently  (e.g.  Monacella 
and  Newman  (1967),  Gadd  (1969),  and  van  Oortmerssen  (1972)).  Of  course  there 
is  more  to  the  solution  of  the  integral  equation  than  just  evaluating  its 
kernel;  however,  direct  numerical  attack  on  this  general  flat-ship  problem 
would  seem  worthwhile,  and  some  effort  is  being  put  into  this. 

The  role  of  the  Kutta,  or  constant- pressure,  condition  at  the  trailing 
edge  is  worth  comment.  There  is  a  degree  of  non-uniqueness  about  the  integral 
equation  (2.11);  the  homogeneous  equation  with  ns=0  has  a  non-trivial  set 
of  solutions.  This  is  illuminated  by  performing  an  indefinite  s-integration 
of  (2.11),  introducing  thereby  an  arbitrary  function  of  x  on  the  left-hand 
side,  say  C(x).  The  resulting  integrated  operator  permits  a  unique  solution, 
the  non-uniqueness  being  now  absorbed  into  C(x).  This  unknown  function  must 


*An  interesting  physical  interpretation  of  this  connection  is  the  statement: 
"Every  planing  surface  is  hydrodynamically  equivalent  to  some  hovercraft." 
The  equivalent  hovercraft  does  not,  of  course,  have  a  uniform  base  pressure. 


somehow  be  determined  by  the  requirement  that  P(x,s)  vanishes  at  the 
trailing  edge.  Physically,  this  indeterminateness  is  equivalent  to  a 
degree  of  indeterminateness  about  the  vertical  location  of  the  hull,  and 
indeed  at  infinite  aspect  ratio  (^  =  0)  ,  C  is  a  constant,  reflecting 
bodily  upward  or  downward  shift  of  the  original  given  foil  relative  to  the 
undisturbed  free  surface  at  infinity. 

The  zero-  and  infinite-Froude-number  limits  of  (2.11)  are  of  interest. 

In  the  zero-Froude-number  case  we  obtain  simply 

P(x,s)  =-pgn(x,s),  (2.12) 

i.e.,  the  appropriate  pressure  is  hydrostatic.  This  is  the  apparent  basis 
for  the  original  flat-ship  formula  of  Hogner  (see  Havelock,  (1932))  which  is 
however  inconsistent,  if  used  in  a  wave-resistance  calculation  at  finite 
Froude  number.  At  infinite  Froude  number,  the  integral  equation  reduces 
exactly  to  that  of  aerodynamic  lifting-surface  theory,  so  that  the  ship  is 
equivalent  to  a  lifting  wing  with  camber  surface  y=n(x,s).  The  role  played 
by  the  Kutta  condition  is  mathematically  the  same;  it  eliminates  a  degree 
of  non-uniqueness  in  the  general  solution  of  the  integral  equation. 

The  analogy  between  the  flat-ship  theory  and  lifting-surface  theory, 
which  becomes  an  exact  equivalence  at  g=0  ,  illustrates  a  disturbing 

feature  of  the  low-aspect-ratio  flat-ship  theory,  namely  that  we  shall  not 
in  general  be  able  to  satisfy  the  Kutta  condition  once  the  low-aspect-ratio 
approximation  has  been  made.  That  is,  the  pressure  predicted  by  the  low- 
aspect-ratio  theory  at  the  edge  of  the  transom  stern  will  not  in  general  be 
atmospheric.  This  would  be  a  most  unfortunate  conclusion,  were  it  not  for 
the  fact  that  low-aspect-ratio  wing  theory  also  suffers  from  this  deficiency, 
yet  nevertheless  has  proved  useful.  What  presumably  happens  is  that  in  a 
small  neighborhood  of  the  trailing  edge  there  is  a  rapid  change  of  pressure 
back  to  atmospheric.  The  hope  is  that  this  occurs  over  a  dynamically- 
insignificant  portion  of  the  total  hull  and  has  no  significant  upstream 
effect.  Some  work  has  been  done  (e.g.  Rogallo,  (1970))  or.  the  corresponding 
aerodynamic  problem. 


3.  Derivation  cf  the  Loti-Aspect- Ratio  Flat-Ship  Integral  Equation 

We  now  assume  that  the  hull  has  a  low  aspect  ratio,  i.e.  that  it  is 
slender,  in  the  sense  that  its  beam  B  is  much  smaller  than  its  length  L, 
say  B=0(£)*  L.  Note  however  that  there  is  a  definite  heirachy  of  smallness 
in  this  problem,  thus 

"Draft«Eeam«Length" . 

The  case  when  the  draft  and  beam  are  comparable  gives  ordinary  slender-ship 
theory,  as  in  Tuck  (1964)  for  low-to-moderate  Froude  numbers,  and  Ogilvie 
(1967)  for  moderate- to«high  Froude  numbers. 

In  the  present  case  we  are  going  to  treat  moderate -to- high  Froude  numbers, 
such  that 

v  =  (3.1) 

U2B 

is  of  order  unity.  This  means  that  the  conventional  length-based  Froude 
number  is  large,  specifically 

f  =  — —  =  0(/l/b)  =0(e“i5).  (3.2) 

t/gL 

This  is  the  regime  treated  by  Ogilvie  (1967)  and  by  Maruo  (1967) . 

In  fact  the  appropriate  integral  equation  can  be  obtained  by  specializing 
Ogilvie 's  (1967)  inner  problem,  for  the  case  of  small  draft/beam  ratio. 

Ogilvie' s  (1967)  general  problem  requires  solution  of  a  non-linear  two- 
dimensional  free-surface  problem  in  each  cross-section.  The  small-draft 
approximation  linearizes  this  problem  and  can  lead  to  the  same  integral 
equation  as  is  obtained  by  the  reverse  procedure,  i.e.  of  "small  a,  then 
small  e,"  rather  than  "small  £  ,  then  small  ot". 

Since  the  body  is  slender,  we  expect  as  usual  to  have  to  solve  a 
two-dimensional  problem  in  the  (x,y)  cross-flow  plane,  i.e.  dropping  $ss 
from  the  Laplace  equation  (2.10)  to  give 

't’xx  +  ^yy  =  ®  (3.3) 

In  the  Froude-number  range  in  which  V=0(1)  it  is  clear  that  both  terms  in 

3  — 

the  free-surface  condition  (2.9)  must  be  retained,  since  -r—  =  O (I-  M  and 

a  ds 

37  =  0(B'1)- 
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If  we  temporarily  define  a  "pseudo-time"  co-ordinate  t  by  the  equation 
s  =  Ut,  (3.4) 

the  free-surface  condition  (2.9)  becomes 


+  4»tt  =  0  (3.5) 

which  is  identical  to  the  usual  unsteady  linearized  free-surface  condition 
for  water  waves.  Thus,  since  <j>  now  satisifies  (3.3),  not  (2.10),  we  can 
use  any  solution  for  unsteady  two-dimensional  linearized  water  waves,  replacing 
t  by  s/U  . 

The  solution  of  most  direct  use  is  again  that  of  a  pressure  distribution 
P(x,s)  over  the  free  surface.  This  is  now  to  be  interpreted  as  "time"- 
varying  pressure  distribution  imposed  on  a  segment  |x|  <  b(s)  of  the 
axis  y-0,  whose  width  2b (s)  also  varies  with  "time". 


The  solution  is  given  by  Wehausen  and  Laitone  (1960,  p.  615) .  It  is 

convenient  to  write  it  in  terms,  not  of  the  velocity  potential  4>(x,y,s), 

but  rather  of  its  conjugate,  the  stream  function  tjj(x,y,s)  .  In  fact,  since 

from  now  on  we  shall  be  concerned  only  with  y=0,  for  brevity  we  write 

<Mx,s)  for  >J>(x,0,s)  •  Thus 

s  b  (O) 


pui{>(x,s) 


d5P(S,a)  K(x-5,s-a), 


(3.6) 


0  -b(0) 


where 


with 

and 


K(x,s)  =  - 


dAsinXx 


=  TTx  F’<“> 


U) 


=  3iif 


4U5|x|  ' 


0) 


F'  (CO)  =  1 


dK 


sin  (C2-co2) 


_d_ 

do d 


f 


d^  cos(^2-oo2) 


(3.7) 


(3.8) 

(3.9) 

(3.10) 
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The  function  F' im)  can  be  expressed  in  terms  of  Fresnel  integrals, 


F’  (u>)  =  1  +  2U)J~  ji 
- 1  *  2“§  1 


COS  ti)  S 


co  so.'2  sinu2 


(Jj  W)  “  sin  w2  C 

sinw2  jJ  .  . , 

2  f,/?a*'3  , 


<Jj  -)  (3-> 


(3.12) 


where  S,C  are  Fresnel  integrals,  and  f(C)  is  an  auxiliary  function 
(Abramowitz  and  Stegun  (1964),  p.  300).  The  function  F'((J)  has  convenient 
series  and  asymptotic  expansions,  respectively 


F'  (aj)  =  1  +  Yl 


(-4W4  )’ 


(3.13) 


m=l  1*3*5*...  (4m-3) • (4m-l) 


which  can  be  used  for  small  in  ,  and 


F'  (o)  =  j^-  (i)(cosw2-sinw2)  -  ^ 


2  r  1  *  3  ‘  5  *  . . .  (4m-3) • (4m-l) 


.  .  4  .  m 

(-4w  ) 


(3.14) 


which  can  be  used  for  large  w  ,  for  a  suitable  stopping  point  m^. 

The  kernel  K  and  function  F'  (u>)  occur  also  in  classical  Cauchy- 
Poisson  problems  (e.g.  Lamb  (1932),  p.  384)  and  the  physical  description  of 
the  spreading  waves  produced  is  well-known.  Indeed  one  can  view  the 
representation  (3.6)  as  resulting  physically  from  a  "time"  history  of 
pressure  pulses,  the  pulse  P(x,s)  at  "time"  a=s  being  applied  in  order 
to  cancel  out  instanteously  the  spreading  waves  produced  at  earlier  "times" 
C<s.  Our  aim  is  to  choose  P(x,s)  so  that  the  stream  function  which  is  left 
over  after  this  cancellation  correctly  satisfies  the  hull  boundary  condition. 
Somewhat  similar  ideas  were  used  by  Cummins  (1956) . 


The  boundary  condition  (2.5)  is  written  in  terms  of  <J> 
using  the  Cauchy-Riemann  equation  <J>  =-i|i  ,  we  have 

y  * 


However, 


iMx,s)  =  -U  j*  ng(S,s)  d £ 


(3.15) 


Note  that  we  have  used  the  natural  antisymmetry  condition  ij»=0  at  x=0. 
Thus  once  the  hull  shape  0(x,s)  is  given,  iMx,s)  may  be  treated  as 
a  known  function,  and  (3.6)  then  represents  an  integral  equation  to 
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determine  the  unknown  pressure  P(x,s).  For  example,  if  the  "ship"  is 
a  flat  plate*  at  an  angle  of  attack  a  then 

n(x,s)  =  -as  ,  (3.16) 

and  we  have  immediately 

*Mx,s)  »  Uas.  (3.17) 

The  analytical  character  of  the  integral  equation  (3.6)  is  of  some 
interest.  The  equation  is  of  Fredholm  character  w'.th  respect  to  the 
(space-like)  variable  x  with  dummy  £  ,  and  of  Volterra  character  with 

respect  to  the  (time-like)  variable  s  with  dummy  a  (Tricomi  1957) . 

This  means  physically  that  information  at  all  values  of  £  is  needed  to 
determine  the  solution  at  any  x,  whereas  only  information  at  further  forward 
stations  o<s  is  needed  to  determine  the  solution  at  a  particular  station 
s.  We  may  hope  to  so) ve  the  equation  in  the  s-dimension  by  a  time-stepping 
or  marching  process,  proceeding  systematically  from  bow  to  stern,  as  in 
an  initial-value  problem  for  a  differential  equation.  But  at  each  station 
s  we  must  expect  to  solve  the  Fredholm  integral  equation  with  respect  to 
x  in  a  manner  more  like  a  boundary -value  problem  for  a  differential 
equation. 

Furthermore,  the  Fredholm  equation  in  the  x  direction  is  singular. 

This  is  most  apparent  at  g=0,  where  u>=0  and  F'(0))  =  F'(0)  =  1.  Thus 
at  g=0  ,  K(x-£,s-0)  =  ^ »  and  the  ^-integration  is  to  be  interpreted 
in  the  sense  of  Cauchy.  In  fact  for  any  g^O  and  O^s  the  singularity  is 
in  a  sense  worse  than  the  simple  Cauchy  pole,  for  as  (V*00  we  have  from 
(3.14) 

F- («)./!«  (cosw2-sinco2) .  (3.18) 


*0r,  in  fact,  any  hull  differing  from  that  given  by  (3.16)  by  addition 
of  a  function  of  x  alone,  for  example,  a  triangular  section  with 
a  constant  deadrise  angle  qualifies.  Only  the  longitudinal  slope  T) 
is  hydrodynamically  significant. 
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Thus  as  Ctx, 


K(x—rJts-o)-* 


,OL 

^8mj2 


(s-a)  ix-c,) 


[g(s-C) ‘  71 
4u2(x-^)+  4 


■V2_0|2is-C) 


(3.19) 


Hence  if  a/s  ,  the  kernel  function  behaves  like  a  "-3/2"  power 
multiplied  by  a  rapidly  oscillating  function,  as  C-*x.  This  behavior 
may  be  expected  to  cause  some  degree  of  numerical  difficulty,  and  does. 


Instead  of  tackling  the  integral  equation  for  the  pressure  P 

itself  directly,  it  is  somewhat  more  convenient  to  work  in  terms  of  a 

function  Q  whose  s  derivative  is  P,  namely 

s 

Q(x,s)  =  J  P(x,Cf)  do  ,  (3.20) 

-.00 


Although  the  lower  limit  of  (3.20)  is  written  as  "-cd",  it  may  equally  well 

be  replaced  by  zero,  or  in  fact  by  sq(x)  ,  where  sq(x)  is  the  station 

s  at  which  x=b(s)  i.e.,  the  function  s  (x)  is  the  mathematical  inverse 

o 

of  the  function  b(s) .  This  is  because  P^O  outside  the  hull  projection  on 
the  plane  y=0. 

The  function  Q(x,s)  is  of  course  the  loading  on  a  unit-width  strip 
of  the  hull  at  offset  x,  extending  from  the  leading  edge  to  station  s. 

Hence,  for  example,  the  total  lift  force  in  the  y  direction  is  obtained 
in  terms  of  the  values  of  Q  at  the  trailing  edge  s=L,  namely 

b  (L) 

Fy  =  f  Q(x,L)  dx  .  (3.21) 

-b(L) 


More-complicated  formulae  involving  Q  at  all  stations  s  apply  to  the 
pitching  moment  and  the  drag.  At  infinite  Froude  number,  Q  is  proportional 
to  the  velocity  potential  <t>;  specifically 

0  =  -PU<!>.  (3.22) 


However,  as  is  clear  from  the  boundary  condition  (2.8),  no  such  identification 
is  possible  if  g/0. 


On  substituting  P(£,0)=Qc(£,a)  in  (3.6),  and  integrating  by  parts  with 
respect  to  o  ,  we  have 


pUiJ;  (x,s) 


b  (0) 

[  d CQ(r,,0)  K (x-£ , s-0) 
-b(0) 


s  b  (a) 

-f  50  J 


o  -b(a) 


d^Q(?,a)Ka(x-C,s-a) 


(3.23) 


b(s)  s  b(C) 

d£&(£,s)K(x-£,0)  -  J  dO  J  dCe(£,0)K0(x-£,s-a)  , 

-b(s)  0  -b(a) 


using  Q=0  at  x=b(s). 


The  first  term  of  (3.23)  involves  K(x-£,0)  which  is  simply  -cy ■ 

Thus  this  term  must  be  interpreted  in  the  sense  of  Cauchy,  and  takes  the 
form  oi  a  finite  Hilbert  transform  (Tricomi,  (1957),  p.  173)  which  we 
write  symbolically  as 


(ait)  ?  J  ■ 


Thus  (3.23)  becomes 


b(s) 

! 

-b(s) 

s  b(O) 

f  -J 

0  -b(0) 


(3.24) 


puifi(x,s)  =^3(s)C(x,s)  -  j*  do  J  dCQ^ojK^x-^s-a) 


(3.25) 


Equation  (3.25)  is  the  principal  integral  equation  we  shall  attempt  to 
solve.  The  kernel  %  may,  after  some  manipulation,  be  written  in  the  form 


Ko(x-S,s-0)= 


2 

ti(s-a) 


3 


F'  (U>) 


(3.26) 


where 

2  _  g(s-a) 2 

w  "  4U 2 1 x-C'I  *  (3.27) 

Equation  (3.25)  agrees  with  the  result  of  integrating  Maruo's  (1967)  equation 
l 57)  with  respect  to  (our)  £  ,  from  £-0  to  £=x.  Table  1  shows  the 

equivalence  of  the  various  symbols  used.  Note  that  the  function  F(x,y) 
used  by  Maruo  in  his  equation  (58)  et  seq.  was  never  defined,  but  is 
related  to  our  K(x,s)  .  Maruo's  equation  (56),  when  similarly  integrated 
with  respect  to  (our)  £  ,  also  agrees  with  our  equation  (3.6). 


TABLE  1.  EQUIVALENCE  OF  SYMBOLS 


THIS  PAPER 

MARUO  (1967) 

S,CJ 

x+£,x'+£ 

x,£ 

y#y' 

Y 

! 

z 

n(x,s) 

f(x,y) 

! 

P(x,s) 

pUY  (x,y) 

Q(x,s) 

pU2g (x,y) 

-n_(x,s) 

5 

w(x,y) 

^(x,s) 

ryto(x,y')dy' 

J0 

b(s) 

bs (x/£) 

B  =  2b  (L) 

2b 

L 

2fc 

F 

(2X0)'l/2 

V 

2X0£/b 

Fy 

L 

4. 


The  High- Froude-Nurrbev  Limit 


The  limit  g-*-0  may  be  carried  out  either  on  (3.6)  or  (3.25).  In 
any  case,  if  g=0  the  kernel  K(x-£,s -a)  =  is  independent  of  s 

and  a.  Hence  in  (3.25) ,  K0=O  and  the  integral  equation  reduces  to 

2$bQ(x,s)  =  pUtMx,s).  (4.1) 

There  is  now  neither  upstream  nor  downstream  influence  of  the  loading  at 
one  station  on  another,  and  the  problem  is  solved  immediately  by  inversion 
of  the  finite  Hilbert  transform,  using  the  inverse  Hilbert  transform 
operator  defined  symbolically  by 

'/fb“1  *  -(b2-x2),/5?t(b2-x2)“,/2  t  (4.2) 

(Tricomi,  1957,  p.  ]79) 

Thus 

Q(x,s)  =  -pu(b2-x2)l/2(^(b2-x2)"l/2^(x,s) 

pU(b2  (s)-x2)  1/2  jp(s)  d^(g,s) _ 

*  -b(s)  (b2(s>-S2)l/2(x-S)  ’  (4.3) 

Normally  the  inverse  operator  dj-  1  is  not  uniquely  defined,  and  to 

1  b 

any  solution  such  as  (4.3)  we  must  add  a  multiple  of  the  function 
(b2-x2)  1 ^2  whose  Hilbert  transform  vanishes  (Tricomi,  (1957),  p.  174). 
However,  we  can  exclude  this  possibility  in  the  present  case,  since  this 
would  generate  velocity  and  pressure  distributions  with  inverse  3/2  power 
singularities  at  the  leading  edge  x=b(s).  In  order  to  retain  only 
integrable  (inverse  square  root)  pressure  singularities,  we  must  require  a 
square-root  zero  in  Q,  leading  to  the  solution  (4.3). 

Since  Q  is  proportional  to  the  velocity  potential  when  g=0  ,  the 
solution  (4.3)  could  also  have  bean  obtained  directly  from  the  boundary- 
value  problem  with  4>=0  as  the  free-surface  condition,  and  simply  expresses 
the  fact  that  conjugate  harmonic  functions  such  as  <j>  and  \p  are  Hilbert 
transforms  of  each  other  on  the  x-axis.  This  solution  is  of  course  well 
known  in  aerodynamic  low-aspect-ratio  wing  theory  (see  e.g.  Newman  &  Wu, 
(1973) .  The  solution  for  a  flat  plate  is  simply 
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Q(x,s)  =  pU2ajb2 (s)-x2  ,  (4.4) 

the  usual  "elliptic  loading  distribution". 

We  are  here  interested  rather  in  the  first  correction  term  to  Q, 
resulting  from  finite-Froude-number  effects.  That  is,  we  seek  an 
asymptotic  expansion  for  small  g  (or  more  correctly  for  small  values  of 
the  appropriate  normalized  gravity  parameter  v=  2—) ,  which  begins  with 

oo  U 

a  term  Q=Q  given  by  (4.3).  Maruo  (1967)  has  performed  such  an  analysis 
on  the  lift  coefficient  in  a  special  case  and  has  proved  that  the  first 
correction  to  the  infinite-Froude-number  lift  is  a  factor  of  order  V  . 

That  is,  the  asymptotic  expansion  at  least  begins  like  a  Taylor  series  with 
respect  to  V  . 

A  logical  procedure  for  constructing  this  expansion  is  by  successive 

approximation,  i.e.  since  (4.3)  resulted  from  dropping  the  last  term  of 

(3.25)  entirely,  the  first  correction  to  (4.3)  is  obtained  by  substitution 
00 

of  Q  into  this  particular  term.  Thus  if  we  put 

Q  =  Q°°  +  Q1  (4.5) 

where  Q^-K)  as  g  or  vK),  we  have 

s  b(0) 

^b(s)Ql  =  |  d0  [  dCQ*(£,o)K0(x-C,s-o).  (4.6) 

0  -b  (0) 

While  (4.6)  may  be  a  useful  formula  as  it  stands,  we  can  simplify  further, 
since  K  itself  still  depends  on  gravity  g.  But  if  for  example  we  were  to 
use  the  (truncated)  series  (3.13)  to  estimate  for  small  g,  we  should 

only  obtain  terms  of  0(g2),  not  0(g)  as  expected  from  Maruo's  analysis. 
Furthermore,  the  resulting  integrals  would  diverge  because  of  a  non-integrable 
singularity  at  £=x  .  The  highly-oscillatory  behaviour  of  the  kernel  near 
£=x  ,  as  indicated  by  (3.19),  suggests  that  the  limit  g-K)  needs  special 

treatment,  and  it  is  clear  by  analogy  with  the  method  of  stationary  phase 
that  only  the  neighborhood  of  £=x  contributes  significantly  to  the  integral 
(4.6),  to  leading  order. 

Thus  we  expand  Q*  in  a  Taylor  series  about  £=x  ,  giving 
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vQ1  =  J  d0  J  <3C|q“(x,o)  +  (£>x)CMx,a)  +...jKo 

0  (nbhd 

of  x)  (4.7) 

However,  £  can  only  take  values  near  to  x  for  so(x)<0<s  where  sQ(x) 

is  as  defined  below  equation  (3.20).  Further,  in  view  of  (3.26),  the 

term  of  (4.7)  in  Q°°(x,o)  integrates  to  zero,  and  we  are  left  with 

s  x+®° 

&bo’  *  -f  [  (4.8) 

Sq  (x)  X-00 

s  00 

=  _£C  da  22si£i£l  f  dxCF*  (a»  -1]  (4.9) 

it  j  s-a  d  » 

sQ(x)  0 

2 

where  X=U“XI  311(3  w2  -  1^2^'"  *  011  chan9ing  the  variable  of 

integration  from  X  to  w  ,  we  have 

^'bQ‘=|  $r  J  do(s-o)Q”(x,o)  f  F’(g~1-  da).  (4.10) 
sQ(x)  0 


The  last  integral  with  respect  to  u)  is  a  pure  constant,  taking  the  value 
"-2tt".  Thus, finally, 

s 

^b(s)Q1(x,s)  =  -  J*  do(s-o)Q”(x,a) .  (4.11) 

sQ(x) 


The  procedure  for  computing  the  leading-order  gravity  effects  on 
the  flow  (  and  in  particular  on  the  pressure  distribution)  is  thus  to 
compute  first  the  infinite-Froude-number  solution  Q°°(x,s)  by  (4.3)  , 
substitute  into  (4.11),  and  then  take  a  further  inverse  Hilbert  transform 
to  find  Q* (x,s)  .  Except  in  very  special  cases  this  procedure  may  be 
nearly  as  difficult  as  direct  numerical  solution  of  the  integral  equation 
(3.25).  However,  we  note  that  the  leading-order  dependence  on  gravity  g 
is  linear,  as  is  Maruo's  (1967)  estimate  for  the  lift  coefficient  of  a 
flat  delta  wing. 
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In  the  special  case  of  the  flat  plate  we  may  proceed  a  little  further. 
Thus  on  use  of  (4.4)  for  Q°°  in  (4.11),  we  have 


s 

^Q1  =  4gpax  £ 


s0(x) 


da(s-a) 
jb*  (0)  -x' 


(4.12) 


It  should  be  observed  that  the  denominator  vanishes  at  the  lower  limit 
0=so(x)  of  the  0  integration.  In  the  further  special  case  of  a  triangular 
waterplane 


b(s)  =  As, 

we  can  integrate  (4.12)  explicitly  to  give 


(4.13) 


#b q1  =  4gp  x  x  s*° g(s+/s2-  fr)  -  j  /s2-  XT  _s*°9  HT 


(4.14) 


Although  no  doubt  the  inverse  Hilbert  transform  could  now  be  obtained 
to  generate  Maruo's  solution,  our  purpose  here  is  rather  to  observe  the 
remarkable  introduction  of  singular  behaviour  along  the  center  line  x=0, 
as  evidenced  by  the  term  of  (4.14)  in  "£og|x|".  This  behaviour  is 
characteristic  of  pointed  flat  plates.  For  example,  if  we  consider  the 
more  general  class  of  waterplanes  whose  behaviour  near  the  bow  is  of  the 


b(s)  =  As11 


(4.15) 


for  some  positive  exponent  n,  then  the  behavior  of  the  integral  (4.12)  as 
x+0  is  of  the  form  of  an  analytic  function  of  x,  plus  a  contribution  of 
the  order  of  s.x1//n  ,  n^l,  or  s.x  iog  x,  n=l. 

Thus  for  all  n>l,  the  slope  of  the  graph  of  the  function 
against  x  is  infinite  at  x=0.  The  function  is  of  course  an  odd 

function  of  x. 

The  corresponding  result  for  Q1  itself  is  that  Q1  behaves  like 
an  analytic  even  function  of  x,  plus  a  contribution  of  the  order  of 
s . | x | 1 for  all  n?l/2,  and  s.xz£og|x|  for  n=l/2.  Specifically  we 


Q1 (x,s)  =  Ql (0,s)  + 


(0  (x2),  0<n<l/2 
s.O  (x2£og|x| ) ,  n=l/2 

s.0(|x|l/n),  n>l/2  . 


(4.16) 
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It  should  be  noted  that  since  P=£s  »  and  the  singular  terms  of  (4.16) 
are  linear  in  s,  (4.16)  indicates  that  the  pressure  distribution  P 
itself  has  the  same  singular  structure,  with  the  strength  of  the 
singularities  invariant  along  the  length  of  the  ship. 

Thus  for  0<n<l/2  ,  the  pressure  is  well  behaved,  for  n=l/2 

(blunt  parabolic  waterplane)  the  lateral  pressure  gradient  vanishes 
at  x=0  but  the  lateral  curvature  of  the  graph  of  pressure  against  x 
is  infinite,  and  similarly  for  l/2<n<l  .  For  n=l  (triangular  water- 

plane)  the  lateral  pressure  gradient  is  discontinuous  but  finite  at 
x=0  while  for  all  n>l  the  lateral  pressure  gradient  is  infinite  at 
x=0.  That  is,  there  is  a  very  sharp  but  finite-magnitude  pressure  peak 
along  the  center  line  of  any  sharply- pointed  flat  plate. 

The  above  result  is  of  course  essentially  a  gravitat icnaj  effect, 
and  stands  in  sharp  contrast  to  the  smooth  behaviour  of  the  elliptic 
loading  (4.4)  in  the  gravity-free  case.  The  singularity  is  presumably 
due  to  the  profound  effect  of  the  diverging  waves  generated  at  the 
extreme  bow,  whose  wave  length  tends  to  zero  along  the  track  of  the 
bow,  irrespective  of  Froude  number  (Ursell,  (I960)).  Although  the  above 
analytic  conclusions  were  obtained  from  a  high-Froude-number  expansion, 
it  is  probable  that  the  essential  character  of  the  singularity  is  the 
same  at  all  Froude  numbers,  and  the  numerical  solution  of  section  7 
tends  to  verify  this.  Experimental  verification  is  eagerly  awaited. 


5.  A  Similarity  Solution 

We  seek  in  the  present  section  a  solution  P(x,s)  of  the  integral 
equation  (3.6)  which  has  the  same  basic  shape  at  all  stations  s.  That 
is,  the  pressure  distribution  at  any  one  station  s  is  obtained  from 
that  at  any  other  station  by  a  simple  scaling  of  P  and  x .  The  x-wise 
scale  is  obviously  b(s),  and  it  is  easy  to  see  that  the  only  possible 
multiplicative  scale  on  P  is  some  power  of  s.  Thus  we  seek  a  solution 
of  the  form 

P(x,s)  =  sYP (fr  *sy  )  ,  (5.1) 

for  some  constant  y  and  some  function  P(x)  of  a  single  variable 
x. 

Of  course  we  have  no  guarantee  that  such  a  solution  ever  exists, 
and  in  particular  could  not  expect  it  to  exist  without  some  special 
restriction  on  b(s).  In  fact, we  shall  now  show  that  (5.1)  is  a  valid 
solution  if  and  only  if  the  waterplane  consists  of  cusped  parabolas, 
i.e.  if 

b(s)  *  Xs2  (5.2) 

for  some  constant  X=B/2L2  .  The  body  shape  n(x,s)  and  thus  the 

stream  function  *Mx,s)  will  also  possess  a  similarity  character,  and 
we  shall  verify  that,  for  example, 


>Mx,s) 


sY‘^(b k)  } 


(5.3) 


The  above  are  the  most  general  assumptions  which  permit  a  similarity 
solution.  To  verify  that  they  are  consistent,  and  to  find  the  equation 
satisfied  by  P(x)  ,  we  substitute  (5.1)  and  (5.3)  into  the  integral 
equation  (5.6). 


Thus  we  have 


s  1 

dC0Yb(0)j*  dCP(C)K 


(x-£b (0) ,s-o) 


(5.4) 


on  substituting  £=£b(0)  .  Putting  x=xb(s)  ,  and  using  the  form 

(3.8)  of  K  gives 
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pUs^+1\Mx)  = 


if  Y  f 

71  \  J  ^ 


(5.5) 


where 


V  =  b(s)/b(o) 


(5.6) 


and 


,2  _  q(s-0)7b(0) 

W  “  4U2 | yx-£ | 


(5.7) 


If  we  put  o=st  in  (5.5),  we  have 

1  1 

<"♦<*>  -  ?  f  «!«>  f  d-fer 

-1  0 


V  (0!) 


(5.8) 


where  now 


and 


V  =  b(s)/b(st) , 


.2  gs2  (1-t)  2/b(st) 

“  =  . 4U^jix-?| - * 


(5.9) 


(5.10) 


So  far,  we  have  not  made  the  assumption  (5.2).  In  order  that  the 
original  assumption  (5.1)  be  valid  it  is  necessary  that  (5.8)  represent 
an  integral  equation  for  P(x)  ,  i.e.  that  it  be  independent  of  the 
station  co-ordinate  s.  The  parameter  y  is  independent  of  s  if  and 
only  if  b  is  proportional  to  some  power  of  s  ,  while  the  parameter  d< 
is  independent  of  s  if  and  only  if  that  power  is  exactly  2.  Thus  if 
(5.2)  is  satisfied,  we  have 

y  =  t~2  (5.11) 

) 


and 

4XU2,x-t2^|  2  |X-t2£|  • 

Thus,  finally,  the  integral  equation  for  P(x)  can  be  written 

1 

PUiMx)  =  f  d^P(C)K(x,Q 

-1 


(5.12) 


(5.13) 
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where  1 

1  P  *4-  4-y+2 

— §)  F’(a,)  ' 
0 

with  u  determined  by  (5.12). 


(5.14) 


The  task  of  solving  the  one-dimensional  integral  equation  (5.13) 
would  appear  to  be  considerably  simpler  than  that  of  solving  the  general 
two-dimensional  equation  (3.6).  However,  the  kernel  K  is  obviously 
extremely  complicated,  and  highly  singular  as  •  Hence  numerical 

similarity  solutions  have  so  far  been  obtained  only  indirectly,  via 
the  qeneral  equation,  and  will  be  presented  in  section  7. 

The  quantity  £(x,s)  whose  s  derivative  is  P  also  obeys  a 
similarity  law,  of  the  form 

Q(x,s)  =  sY+1  Q(^y  >  •  (5.15) 


On  differentiation  with  respect  to  s,  we  establish  the  connection 


P(x)  =  (Y+l)  Q(x)  -  2x  Q‘ (x)  (5.16) 

•s,  <v  <v  “V 

between  the  similarity  profiles  of  P  and  Q  .  An  integral  equation  for 

Q  may  also  be  obtained  by  substitution  of  (5.15)  into  (3.25), 

namely 

1 

pUijJ(x)  «#jQ(x)  +  j*  dCQ(0-^K1(x,0  ,  (5.17) 

-1 


where 


Kx(x,p 


2  r  dt  tY+1 

,  J  1-t 

o 


[F 1  (u>)  -1)  , 


(5.18) 


0)  being  given  by  (5.12)  again. 

The  allowed  shapes  of  the  hull  are  of  interest.  Clearly  we  are 
allowed  only  the  very  specific  cusped  waterline  prescribed  by  (5.2)  . 
However,  considerably  greater  latitude  is  allowed  in  the  shape  of  the 
cross-sections  and  some  latitude  is  allowed  in  the  longitudinal  profile. 
Thus  it  is  clear  that  the  hull  function  n  also  has  a  similarity 
character,  with 
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n(x,s) 


sYn( 


(5.19) 


The  arbitrary  exponent  Y  therefore  characterizes  the  longitudinal 
profile  or  keel  line,  which  is  straight  if  Y=1  »  blunt  if  Y<1  and 

cusped  if  y>l*.  The  case  of  a  flat  plate  is  included,  with  Y=1 
and  n<x)=  -a=constant.  More  generally,  any  cross  section  shape  defined 
by  n<x)  is  al.’o'-sl,  but  the  hull  shape  must  of  course  be  similar  for 
all  stations,  according  to  (5.19).  The  connection  between  the  shape 
function  n(x)  and  stream  function  if>(x)  may  be  obtained  from  (3.15), 
and  we  have 


x 


ijj(x)  =  2XUxn(x)  -  Xu(Y+2) 


s 


n(Od§  . 


(5.20) 


A  physical  justification  of  this  similarity  solution  may  be  attempted 
as  follows.  At  these  high  Froude  numbers  we  are  concerned  only  with 
the  diverging  part  of  the  ship  wave  pattern  near  the  ship's  track, 
since  the  transverse  wavelenqth  27TU2/g  far  exceeds  the  ship  length. 

The  diverging  waves  are  (Ur sell,  (I960))  short  in  wavelength  even  for 
vanishing  gravity,  and  in  fact  their  crests  asymptote  to  the  axis  x=0 
accorling  to  the  parabolic  law,  x~s2 .  Thus  the  growth  of  the  waterplane 
(5.2)  precisely  matches  the  spreading  of  the  diverging  waves. 

Given  this  physical  picture  we  may  speculate  on  the  character  of 
the  solution,  especially  near  the  leading  edges,  for  other  waterplanes. 

For  example,  if  the  waterplane  is  more  highly  cusped  than  (5.2),  i.e. 
n<2  in  (4.15),  the  rate  of  spreading  of  waves  exceeds  the  rate  of 
growth  of  waterplane,  and  at  some  station  the  diverging  waves  must  emerge 
from  beneath  the  hull,  cnanging  fundamentally  the  character  of  the  leading 
edge  singularity  and  hence  the  spray  sheet. 


*  Remarkably,  the  special  case  y=2  allows  similarity  solutions  with  the 
non-linear  free  surface  condition  (2.6). 


6.  Numerical  Procedure 


In  this  section  we  discuss  a  procedure  for  numerical  solution  of 
the  integral  equation  (3.25).  The  program  used  is  quite  unsophisticated, 
and  further  work  is  needed  to  develop  more  efficient  programs.  However, 
the  accuracy  attainable  with  the  present  method  is  satisfactory  for 
some  purposes. 


The  only  numerical  difficulty  in  solving  (3.25)  is  with  the  double¬ 
integral  term.  Routines  for  efficiently  inverting  finite  Hilbert 
transforms  are  easy  to  construct,  so  that  the  first  term  on  the  right  of 
(3.25)  gives  no  trouble.  Notice  that  the  double-integral  term  contains 
all  "time"-history  effects;  that  is,  it  and  only  it  introduces  an 
influence  of  previous  stations  a<s  on  the  pressure  at  the  current 
station  s.  In  this  connection  it  is  important  to  note  that  the  kernel 

K  vanishes  at  the  current  station,  ie.  when  0=s . 
o  — — — — 

We  make  use  of  this  property  in  a  "time"-stepping  procedure,  by 
first  using  the  ordinary  trapezoidal  rule  on  the  0- integration.  Having 
chosen  a  station  spacing  As, we  write 

Qn(x)  =  Q(x,nAs)  (6.1) 

^(x)  =  4>(x,nAs)  (6.2) 

b  =  b(nAs)  (6.3) 

n 


and  approximate 

pity  (x)  =«| k  Q  (x) 
n  '  Dp  n 


n-1  bk 

As  Z  dCQk(C)Ko(x-C,  (n-k)  As) 
k=l  -bk 

(6.4) 


Equation  (6.4)  can  be  written  in  the  form  of  a  recursive  algorithm 
for  Qr,  namely 

2n(x)  ^b1  Rn(x)  (6.5) 

n 


where 


n-1 

R  (x)  =  pUiP  (x)  +  As  5" 
n  n  «-• 


k=l 


J*  d?Qk(C)Ka(x-C,  (n-k)As)  . 

'bk  (6.6) 
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Note  that  R^  is  determined  from  the  known  quantity  ^(x)  »  together 

with  all  Q^(x),  k=l,2, . . .n-1,  which  are  known  at  the  n'th  step. 

The  next  problem  is  evaluation  of  the  ^-integral  in  (6.6). 

We  use  a  very  crude  estimate,  in  which  Q  (fj)  is  taken  as  a  constant, 

n 

say  Q..  ,  on  each  of  2M  segments,  j=±l,±2, . . .±M  ,  where  the  j'th 

segment  is  defined  by 


Ej-i  <  E  ‘  ? ‘  bk  sin  m  •  <e 

Note  that  the  same  number  2M  of  x-wise  segments  is  used  at  every 
station  s,  the  segment  size  increasing  with  waterplane  width  2b^. 

Vsing  the  formula  (3.26)  for  and  integrating  explicitly  with 

respect  to  £  ,  we  have 


(6.7) 


F  (x)  =  QUip  (x)  -  ;■  r 
n  n  TT  *—• 


z  grH? 


(6.8) 


^'j-1 


2  _  q  (As)  2  (n-k)  2 
w  4U2  |X-£| 

We  now  evaluate  (6.8)  at  a  point  x=x^  ,  which  is  approximately  the 
mid-point  of  the  i'th  segment,  namely 

.  .  (i— *s)  tt 

xi  =  bk  sin  2M  ' 


(6.9) 


(6.10) 


obtaining  a  set  of  values  R^n=Rn(xJ  *  Finally,  a  numerical  inverse 
Hilbert  transform,  essentially  evaluating  an  expression  like  (4.3)  by 
the  mid-point  rule  (after  removing  the  Cauchy  singularity) ,  provides  a 
corresponding  set  of  values  of  =  Qn(xJ  •  We  now  proceed  to  station 

n+1  ,  etc.  Note  that  no  matrix  manipulation  (especially  no  inversion) 

is  ever  required  in  this  method. 

Difficulties  arise  because  of  the  highly-singular  nature  of  the 
kernel  near  £=x  ,  as  indicated  by  (3.19).  Of  course  we  never  evaluate 

exactly  at  this  point,  and  in  obtaining  (6.8)  from  (6.6)  have  integrated 
analytically  through  this  singularity.  Nevertheless,  there  is  bound  to 
be  trouble  in  (6.8)  due  to  large  values  of  U)  whenever  the  point  of 
evaluation  x=x^  at  some  station  s  is  close  to  an  end  point 
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of  any  segment  at  a  previous  station  a  .  Physically,  each  end  point 
of  a  segment  looks  like  an  isolated  singularity,  which  leaves  its  own 
"trail"  in  the  form  of  wildly-oscillating  waves  (Ursell,  (I960)).  A 
better  numerical  method  could  be  one  in  which  the  step-function 
character  of  with  £  was  replaced  by  a  smoother  variation, 

thereby  moderating  the  apparent  singularity. 

This  problem  manifests  itself  in  the  form  of  apparently-random 

small  fluctuations  of  R^fx)  as  a  function  of  x,  superposed  upon  a 

"believable"  smooth  wave.  It  is  cured  in  a  not-altogether-satisfactory 

manner  by  two  separate  smoothing  procedures.  In  the  first  place,  we 

test  each  end  point  while  evaluating  the  sum  (6.8)  to  ensure  it  is 

not  too  close  to  the  current  evaluation  point  .  If  it  is  as  close 

to  x^  as  20%  of  the  ith  segment  size,  we  shift  x^  (staying  within  the 

ith  segment)  by  that  20%  amount.  Secondly,  after  complete  evaluation  of 

the  R  (x, )  ,  we  smooth  by  replacing  R  (x.)  by  the  mean  of  R  (x.  ,) 

ni  n  l  n  l-l 

and  V*ltl>. 

The  particular  trigonometric  lateral  spacing  (6.7)  was  chosen  to 
provide  a  sufficient  density  of  segments  near  the  edges  to  counter  the 
rapid  (square-root)  drop  to  zero  of  Q.  In  fact,  explicit  use  is  made  of 
the  nature  of  this  spacing  to  make  the  inverse  Hilbert  transform  most 
efficient,  and  for  example  the  program  reproduces  exactly  the  result 
(4.4)  for  flat  plates  at  infinite  Froude  number.  However,  this  decision 
was  made  before  the  singular  character  of  Q  near  the  center  plane  was 
discovered.  Actually  the  investigation  of  Section  4  was  only  carried 
out  as  a  result  of  the  appearance  of  the  numerical  results,  and  in 
retrospect  it  would  appear  that  a  greater  density  of  points  near  the 
center  plane  would  have  been  desirable. 


7.  Discussion  of  Computed  Results 

Figure  2  shows  results  for  Q/pU2ab(s)  plotted  against  x=x/b(s) 
at  various  stations  s,  for  the  case  of  a  flat  plate  with  the  cusped 
parabolic  waterplane  (5.2).  This  is  the  case  in  which  a  similarity 
solution  exists,  such  that  the  quantity  plotted  should  be  independent 
of  s.  The  results  shown  are  for  M=20  and  a  maximum  value  of  n=20  , 

with  the  speed  chosen  so  v=1.25  .  For  example,  with  a  length /beam 
ratio  of  5.0,  this  would  correspond  to  a  conventional  Froude  number  F=2.0  . 

We  observe  that  at  this  fairly-high  Froude  number,  a  similarity 
profile  is  reasonably  well  achieved  by  about  the  mid-section  of  the 
ship.  In  fact  departure  from  similarity  very  near  the  bow  is  inevitable, 
since  the  program  starts  with  R^(x)  =0  in  (6.5).  That  is,  at  the  very  first 
station  n=l  there  is  an  apparent  infinite-Froude-number  or  zero-gravity 
solution,  irrespective  of  the  actual  Froude  number.  This  is  shown  as 
the  s=0  curve  in  Figure  2,  and  is  simply  the  elliptic  loading 
(4.4).  The  behavior  for  the  first  few  stations  is  quite  erratic,  but 
the  oscillations  apparently  die  out  as  s  increases. 

If  we  now  vary  v  ,  i.e.  vary  the  Froude  number,  we  obtain  the 
family  of  similarity  profiles  shown  in  Figure  3.  These  are  essentially 
plots  of  Q (x)  ,  as  in  (5.15)  with  y=l.  However,  they  are  actually 
obtained  as  in  Figure  2  from  the  general  program  at  s/L=1.0 
Similarity  is  harder  to  achieve  numerically  as  v  increases,  i.e. 
as  the  effect  of  gravity  increases,  especially  near  the  center  plane. 

The  curves  are  dashed  wherever  an  uncertainty  of  more  than  about  5% 
exists,  and  discontinued  altogether  as  soon  as  the  uncertainty  reaches 
10%. 

Corresponding  curves  of  the  actual  pressure  P(x)  could  be 
obtained  from  (5.16),  but  the  necessary  numerical  differentiation  would 
reduce  the  accuracy  of  the  results  unacceptably.  However,  it  is 
clear  that  the  general  character  is  of  a  sharp  but  finite  pressure  peak 
at  the  center  plane,  with  a  pressure  minimum  about  half-way  out,  followed 
by  an  infinite  positive-amplitude  (inverse-square-root)  peak  at  the 
edge  x=l.  This  infinity  corresponds  to  the  leading-edge  spray  sheet. 
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Figure  2:  Similarity  check  for  cusped  waterplane.  Scaled 
distribution  at  various  stations,  for  fixed  v  • 
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Figure  3:  Trailing  edge  loading  distribution  for  cusped  v/aterplane 
at  various  values  of  v  =  gL2/U2B 
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Figure  4  shows  the  loading  Q  at  the  trailing  edge  s/L=1.0  for 
the  case  of  a  triangular  waterplane,  i.e.  for  a  delta  wing,  at  various 
values  of  V  The  results  are  analogous  to  those  of  Figure  3,  but 
the  profiles  are  no  longer  self-similar  with  respect  to  s.  On  the  contrary, 
Figure  4  has  an  alternative  interpretation  as  a  plot  of  scaled  loadings 
at  various  stations  s  for  a  fixed  value  of  v  .  For  example,  at 
v=2.5  the  loading  at  the  mid-section  s/L=0.5  is  precisely  half  of  the 
result  shown  in  Figure  3  for  the  trailing-edge  loading  at  v=1.25. 

In  both  Figures  3  and  4  the  center-plane  singularity  suggested  by 
the  analysis  of  Section  4  is  qualitatively  evident.  Unfortunately, 
program  accuracy  in  this  region  is,  not  surprisingly,  least  satisfactory, 
so  that  we  are  not  able  to  verify  the  differences  in  tne  actual  degree 
of  singularity,  as  indicated  by  (4.16). 

Perhaps  a  more  significant  difference  between  Figures  3  and  4  is 
in  the  strength  of  the  pressure  singularity  at  the  edge  x=l,  which 
appears  approximately  invariant  with  v  for  the  similarity  profiles  of 
Figure  3.  In  the  case  of  the  delta  wing,  however,  there  appears  to 
be  a  real  weakening  of  the  pressure  singularity  as  V  increases,  or 
as  we  move  from  bow  to  stern  at  fixed  v  .  In  fact  for  all  v>2.1 
the  computer  program  predicts  small  negative  loadings  very  near  to  x=l. 
Since  this  implies  a  negative  infinity  in  the  pressure,  it  is  not  a 
physical ly-acceptable  result.  Unfortunately  it  is  hard  with  the  present 
crude  program  to  tell  whether  these  are  genuine  theoretical  predictions, 
or  numerical  errors.  However,  the  fact  remains  that  no  such  negative 
values  are  ever  obtained  in  the  similarity  case  of  Figure  3. 

This  effect  is  anticipated  by  the  discussion  at  the  end  of  Section 
5,  and  we  can  illustrate  it  more  strongly  by  use  of  the  even-blunter 
waterplane  n=l/2  in  (4.15),  ie.  one  which  is  parabolic  in  s  against  x. 
Figure  5  shows  loadings  for  this  case.  There  is  now  no  doubt  from  the 
computer  output  that  for  V>1.1  the  predicted  edge  loading  is  negative. 

What  actually  happens  here  is  not  clear;  what  is  clear,  however,  is  that 
the  present  theory  is  no  longer  valid.  Note  also  from  Figure  5  that  for 
this  blunt  body,  the  center-plane  singularity  has  almost  disappeared  and 
the  pressure  gradient  now  appears  to  vanish  at  x=0  ,  as  predicted  by 

Section  4. 


Figure  4:  Trailing  edge  loading  distribution  for  triangular 
(delta  wing)  waterplane,  at  various  values  of  v  = 
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5:  Trailing  edge  loading  distribution  for  blunt  vaterplane 
at  various  values  of  v  =  gL2/U2B 
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Figure  6  shows  the  variation  with  v  of  the  lift  force  F  , 
computed  according  to  (3.21)  and  scaled  with  respect  to  the  infinite- 
Froude-number  (i.e.  v=0)  value 

F°°  =  I  Pu2a(b(L))2  .  (7.1) 

y  2 

Results  for  all  three  waterplanes  discussed  above  are  shown.  For  the 
triangular  case  only,  comparison  may  be  made  with  Maruo's  (1967) 
very-high-Froude-number  approximation 

=  1  +  0.211V  ,  (7.2) 

clearly  gives  the  correct  asymptotic  behavior  for 

small  v  . 


Vry 
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Figure  6:  Lift  of  various  flat  plates,  scaled  with  respect  to 

the  zero-gravity  limit,  and  plotted  against  v=aL2/U2B. 
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APPENDIX:  Computer  Program  Listing 


MAIN  PROGRAM 


C  LCW  ASPECT  RATIO  FLAT  SUP  THEORY 

C  COMPLIES  ICAUHG  J  (*STA1  ILNkISE  PRESSURE  IMtGRAL)  AND  THE 
C  RESULTING  LIFT  CUEFFIC IEN  T  FUR  SHIPS  OF  SMALL  DRAFT  ANO  BEAM 
C  REFERENCE:  MARJU  1S67,  TUCK  1S/4 
C 

DIMENSION  R( 50 ) ,GQ ( 50  ) » PS  I (  50, 50)  *  0X150 

CCMMCN  XM( 50) *  SXC50) ,M, UPHl OP, A( 50 ) , TT ( 50 ) , G*XX ( 5U ) , Q< 50, 50) * TPOT 

M  =  NUMEER  OF  LATERAL  PGIMS  (OFFStTS , dUT TOCKS)  IN  A  HALF-HIOTH 
OF  THE  SHIP.  NOTE:  SAME  FCR  ALL  STATIONS.  SPACING  IS  SQUARE  ROUT 
dIAiSEC  TUrfAKD  EDGE*  SCALtO  *KT  HALF  SEAM  A(J)  Af  STAflON  J  . 

N  =  NUMBER  OF  STATIONS,  EQUALLY  SPACED 

RE AC  (5,9)  M»N 
9  FORMAT (2110) 

WRI TE (6, 21 )  M,N 

21  FCRMAT(4H1  M=,I4,3H  N*, 14,///, 32H  J  TT(J)  A(J)  PSKJ.K)) 
fcM  *  M 

OPFIOP  =  C.5/EM 
OPHI  =  3.1416  *  DPHIOP 
PCPHI  *  C.5  *  UPHl 
JQ  2  K  =  1,M 
PHI  =  K*CPHI 
XX ( K )  =  SINiPHI ) 

PHI  =  PHI  -  HOPHI 
XM ( K )  *  SINIPHI ) 

S  X ( K )  =  CCSIPHI ) 

DX( K )  =  SX(K)  *  OPHI 
2  CONTINUE 
C 

C  PRESENT  PRUGRAM  GENERATES  HULL  OATA  INTERNALLY  FOR  FLAT  OELTA 
C  WING, HALF  APEX  '* ANGLE’*  =0.1,  ANGLE  JF  ATTACK  ”SLOPE"*U.  I,  SO  THA T 
C  THE  HALF  RATERPLANE  wlUTF  IS  A(J)  =  0.1  *  T  T  ( J)  .  TT  IS  T  HF 
C  STATICN  CGQRU,  GOES  FROM  C  TU  1  .  GENEK At  I  SAT  I CN  TO  MORF  GENERAL 
C  RATERPLANES  IS  EASY;  JUST  REPLACE  DEFINING  STATEMENT  FOR  A(J). 

C  GENERALISATION  TO  OTHER  THAN  FLAT  PLATE  REQUIRES  MURE  EFFORT,  SEE 
C  TUCK  1974,  TU  SET  UP  MATRIX  OF  STREAM  FUNLTICN  PSIlJ»K) 

C 

SLCPG  =  0.1 
ANGLE  =  0.1 
OT  =  i./N 

TPOT  =  0.63661  *  DT 
AREAwP  =  0. 

DO  1C  5  J  =  1,N 
TT ( J )  =  CT  *  J 
A  ( J )  =  ANGLE  *  TT ( J  ) 

ARE ArfP  =  AREAWP  ♦  A(J) 

DO  100  K  =  1  ,M 

100  PSKJ.K)  =  SLCPE  *  XM(K)*A(J) 

105  WRI TE(6, 106)  J , T T( J)  ,  A ( J)  ,( PS  I ( J , K) , K  *  1,M) 

106  FORMAT (16, 12F9.5) 

AREA  i,P  =  2.*0T*(AREAfcP-0.5*A(N)  ) 

WPCC6F  =  0.5  *  AREAWP  /  A(N) 

ASPECT  =  (4.  *  A (N ) **2  )  /  ARE  A*|P 
i3L  =  2.  *  A(N) 
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Main  Program  (continued) 


WRITE(6,61)  ttl. ASPECT, nPCUEF 

61  F0RPAT114H  OtAP/LENGTH  =,F8.4,15H  AsPELT  RAIIO  =,F8.4,19H  V.ATERPLA 
1NE  CCEFF  =,F8.4) 

INPUT  VALUES  OF  G=l/F**2  «  F=F80UDE  NC.  BASEC  ON  LENGTH. 

NEGATIVE  G  HAS  EFFECT  OF  ZERO  G  (INFINITE  F  } 

ZERO  INPUT  G  VALUE  (THAT  IS. BLANK  CARO)  STCPS  PROGRAM 

10  REAU2.ll)  G 

11  FCRPAT  (FiO.S) 

IF(G )  12.13,14 

12  *RITE(6,IS) 

15  FCRPAT  ( 5H1  G*0) 

i»  s  0« 

GC  TC  17 

14  FROLOE  -  1 ./SORT (G) 

«RITF(6,16)  G, FROLOE 
lb  FORMAT (4H1  G-  ,  Fb.3 , 7FFRUUUE- , F  6. 3) 

17  aR I T  £ (6,63 ) 

63  FORMAT  (26H  UU.l)  C(J,2)  ETC  ) 

00  131  J  *  l.N 
CALL  FLAT ( J  ,R) 

CALL  SMOCTH(R) 

OG  1C2  K  =  1 »M 

102  «(K )  *  PSI(J.K)  -  K ( K  ) 

CALL  HIL B 1N( Oti* R ) 

■JC  44  K  *  l.M 
44  C(J  ,K)  *  00  ( K ) 

101  WRITEI6, 103)  (U ( J  ,K) »K-1 ,M) 

103  FORMAT ( 10F10.5) 

CLIFT  *  0. 

00  60  K  *  l.M 

60  CLIFT  *  CLIFT  ♦  Q(N.K)  ♦  OX(K)  *  A(N ) 

CLIFT  *  4.  *  CLIFT  /  AREAtmP 
HR  I TE (6. 62 )  CLIFT 
£2  FORMAT (F 10.5) 

GO  TO  10 

13  STCP 
cND 


End  of  Main  Program.  Listings  follow  for  Subroutines  FLAT,  SMOOTH, 
HILBIN  and  FPRIM2. 
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SUfiROUTlNE  FUHM.KI 
L  EVALUATES  R(J)  ,  AS  IN  TUCK  157* 

DIMENSION  R(50) 

CLMMCN  XMI50)  ,  SXT5UI  ,M,DPhIOP,A(  50),m5C),G,XXl5G),Q(  50, 50), T PD 1 
DU  5  KK  *  I ,  M 
5  R(KK)  =  C. 

IF(M-I)  lOU,  100,101 

101  IMG)  100,100,102 

102  N  =  M  -  I 

XCLO  =  -  A (Nil  *  XM(l) 

UU  3  KK  -  1 ,  M 

X  *  AIM!  *  XMiKKi 
C 

C  THIS  IS  A  CRITICAL  DECISICN  lard.  JXLKir  IS  TEE  critical 
c  closeness  between  the  current  field  point  anc  the  track  cf  a 

C  PRE  VIOL  S  PRESSURE  PCINT.  THE  NUMBER  "0.2"  USED  IS  RATHER  ARBITRARY 
C  VALUES  OF  0.1  AND  C.3  GIVE  SIMILAR  RESULTS. 

UXCRIT  =  J.2  *  (X  -  XCLD) 

SUMCUT  =  0. 

CC  2  J  =  i,N 
T  =  TTINI  t  -  THJi 
SUM  IN  =  0. 

wNUM  =  0.25  *  G  *  T**2 
FKOLD  =  0. 

DC  I  K  -  I »M 

XI  *  A(J)  *  XX  ( K ) 

AX  *  ABSIX-XI ) 

C  HERc  IS  WHERE  JXCR  IT  IS  USED.  IF  "AXMI*  BELOW  OXCRIT,  WE  SIMPLY 
C  REPLACE  IT  BY  UXCKIT,  THUS  CALLING  ThE  FRtSNEL  INTEGRAL  ROUTINE 
C  WITH  A  SUBSTANTIALLY  REOLCEC  ARGUMENT.  A  VERY  RCUuH  TRICK 
IFIAX.LT. OXCRIT)  AX  «  DXCRIT 
WMNLS2  *  WNUM/AX 
wPLLS2  *  WNUM/(X«-XI) 

FKERN  =  FPR IM2 ( WMNUS2  )  -  F PR  I  M2 (WPLUS2 ) 

SUM  I N  *  SUM IN  ♦  0 ( J , K  )  *  ( FKEKN  -  FKOL  U ) 

1  FKCLJ  =  FKEKN 

2  SUMCUT  =  SUMUUT  F  SUMIN/T 
XCLC  =  X 

3  R(KK)=TPOT*SUMQCT 
100  KETIRN 

END 


SUBROUTINE  SMUUTH(K) 

C  SMOOTHS  R  BY  REPLACING  CLD  H  WITH  THt  AVERAGE  UF  ITSELF  AND 
C  *HA  T  WE  GET  BY  LINEAR  INTERPOLATING  BETWEEN  THE  2  NEAREST  VALUES. 

COMMON  X  M ( 50 ) ,  S X ( 50 ) ,M,DPH IOP, A{ 5 J) ,T T ( 5C ) , G ,XX ( 50  ) , G( 50, 50) , TPDT 
DIMENSION  R (50 )  ,  RNI 50) 

Ml  «  M-L 
DO  I  J  *  2, Ml 
JM  »  J  -  1 

JP  “  J  ♦  1 

0*1./  (XM(JP)  -  XM(JM)  ) 

AM  *  (  XM(J)  -  XM(JM)  )  *  0 

AP  *  (  XN(J)  -  XM(JP)  )  *  D 

1  RN(J)  *  AM  *  R(JP)  -  AP  *  MUM) 

RN(M)  «  ((I.  -  XM ( M )  )/(!.-  XM(Ml)  )  )  *  R(Mi) 

00  2  J  *  2 ,  M 

2  K ( J )  *  0.5  *  (  R ( J )  F  PN(J)  ) 

RETURN 

END 


II 
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SUBRQUTINE  HlLBIN  (<J,K) 

C  INVERSE  HILBERT  TRANSFORMER,  USES  MIO  PT  RULE 
C  OUTPUT  Q  AT  SAME  (COSINE)  SPACINGS  AS  INPUI  R 
C  SPACINGS  MUST  BE  COSINE,  THAT  IS  SUUARc  RC2T  8IAS  TO  ENDS 
C  ASSUMES  ANTISYMMETRY  OF  INPUT  R  ,  SYMMETRY  OF  OUTPUT  C  ,  ABOUT 
C  CENTERLINE  X  *  0.  USES  CNLY  POSITIVE  A*S,  BLILOS.  IN  SYMMETRIES. 
COMMON  XM(50),  SX(50),  M,DPHICP 
DIMENSION  H  50)  ,  R(SO),  0(50) 

OC  1  L  -  I ,  M 
SUM  »  0. 

DC  2  K  -  It M 
IF  (  K  -  L  )  3,2,3 

3  FIX)  *  (  ft(L)-aU)) /(XM(L)-XM(K) )  *  (R(L)+R(K))/(XM(L)  ♦XM(  K  ) ) 

2  CONTINUE 

IF(L-l)  4,4,5 

4  F(l)  *  2  «  *F ( 2  )  -  F  ( 3) 

GC  TU  8 

5  JF(L-M)  6,7,7 

7  l-(MJ  *  2.*F(M-1 )  -  r  ( M-2 ) 

GC  TO  8 

6  F ( L )  =  0.5  *  (F(L-l)  ♦  MLOI  ) 

8  OC  9  K  -  I, M 

9  SUM  =  SUM  ♦  F  ( K ) 

1  y(L)  =  OPHIUP  *  SX(L)  *  SUM 
RETURN 
END 


FUNCTION  FPR IM2 (MW  ) 

C  FRESNEL  INTEGRAL  PCUTINE 

DOUBLE  PRECISION  ZZ ,UC , Tfc R , FPR IN 
IF(NW-l6.)  1,2,2 

1  ZZ  =  -  4.  *  RW**2 
FPR I M  *  1. 

TER  *  1. 

DO  3  M  -  1,500 
MM  =  4  *  M  -  1 
00-  MM  *  ( MM- 2 ) 

TER  =  TER  *  ZZ/  00 

FPR I M  =  FPRIM  ♦  TER 

IF(CA8S( TER)  -  0.000001)  4,3,3 

3  CCNTIMJE 

4  FPRIM2  =  FPRIM 
RETUFN 


2 


5 

6 


a  *  SCRT(WW) 

L  =-4.  *  RW**2 
SUM  *  l. 

TERM  =  1. 

DO  5  M  =  1,20 
MM  =  4  *  M  -  1 
0  =  MM  *  ( MM-2 ) 

TERM  *  TERM  *  U  /  Z 

SUM  *  SUM  ♦  TERM 

IF (ABS(TERM)  -  G.OOOOC1)  6,5,5 

CCNT  INUE 

FPR I M^  =  i.  -  SUM  *  1.253314  *  W  *  (COS(NU) 


RETURN 

END 


Reproduced  from 
best  available 


SIN(HM)  ) 


